Landau critical velocity in weakly interacting Bose gases 
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The flow of a uniform Bose gas at speeds greater than the Landau critical velocity, v c , does not 
necessarily destroy superfluidity, but rather need only lead to a decrease of the superfluid mass 
density, p s . Analyzing a weakly interacting Bose gas with a finite range interparticle interaction 
that leads to a Landau critical velocity at non-zero quasiparticle momentum, we explicitly construct 
the (non-uniform) condensate for fluid flow faster than v c and calculate the accompanying decrease 
in p s . We briefly comment on the relation of the physics to other problems in superfluids, e.g., 
solitons, and vortices in Bose-Einstein condensates, and critical currents in superconductors. 
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I. INTRODUCTION 

When a superfluid flows sufficiently rapidly that the 
excitation spectrum at non-zero momentum in the lab 
frame becomes gapless, then according to the Landau 
criterion, superfluidity ceases 0,0. However, in many 
situations, e.g., dilute solutions of 3 Hc in superfluid 4 Hc, 
as well as superfluid 4 He at non-zero temperature, the 
excitation spectrum is indeed gapless. Yet these systems 
remain good superfluids, only with the superfluid mass 
density p s less than the total mass density p. As a conse- 
quence such systems have a non-vanishing normal mass 
density, p n = p — p s as well, and thus exhibit two-fluid 
behavior. On the other hand, a system with a gap in the 
excitation spectrum is not guaranteed to be superfluid, 
as one sees for example in the Mott insulating phase of 
a Bose gas in an optical lattice Q. In this paper we 
study in detail the state of superfluidity in a weakly in- 
teracting Bose gas flowing at a velocity greater than the 
Landau critical velocity. Pitaevskii has laid out the basic 
physics of this situation in the context of rotons in su- 
perfluid 4 He [J; once the rapidly flowing system begins 
spontaneously to create excitations, interactions between 
the excitations -when repulsive - modify the excitation 
spectrum so that the system remains just critical. 

We address here flow of the system past walls at speeds 
v greater than the Landau velocity, v c , The situation in 
which an object moves faster than v c through a super- 
fluid at rest is different: there the object will generate 
excitations as it moves and will experience dissipation, as 
observed in cold atom experiments [f|. The excitations 
propagate away, and do not affect the energy required for 
the object to make a further excitations. 

Explicitly, in a superfluid with an elementary excita- 
tion spectrum e(p) at momentum p flowing with velocity 
v (in the positive x direction) with respect to the sys- 
tem walls, the excitation energy in the frame at rest is 
e(p) + vp x . The Landau criterion states that the super- 
fluid will become energetically unstable when e(p)/p = v 



(for p in the negative x direction, the direction oppo- 
site to the flow). At the minimum v = v c satisfying this 
condition, the group velocity of the excitation, de{p)/dp, 
equals its phase velocity e(p)/p. We consider here a 
weakly interacting Bose gas at zero temperature initially 
fully condensed into a state of flow in the x direction with 
condensate momentum q > 0, and a corresponding order 
parameter in the lab frame, 



tp(x) = \/ne l 



(1) 



with v = q/m. With a zero range interaction, the exci- 
tation spectrum curves upward, and the group velocity 
always exceeds the phase velocity at non-zero momen- 
tum. Thus to study the Landau criterion, we assume 
a finite range interaction g(\r — r '\) between particles, 




FIG. 1: The assumed quasiparticle excitation spectrum as a 
function of momentum. The slope of the dashed line is the 
Landau critical velocity, v c - 
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for which the Bogoliubov excitations of the system with 
momentum p have energy 



e(p) 



ng(p)p 2 



P 

2m 



-,1/2 



(2) 



where n is the density and g(p) is the Fourier transform 
of g(r). We take g(p) positive and decreasing with p. 
Then the phase and group velocities of an excitation are 
equal when 



nm dg{p) d(gn) 



p dp d(p 2 /2m) 



(3) 



We assume that this condition is satisfied at the criti- 
cal momentum p c ; thus the flow velocity at the point at 
which the Landau criterion is satisfied, is 



t-'c = [{Pcf2m) 2 + ngi/m 



1/2 



(4) 



where g\ = g(p c ), and more generally for integer v, g v = 
g{vp c ). The condition that the velocity, q/m, exceeds v c 
can be written as 



q 2 > p 2 /4 + mng\. 



(5) 



The excitation spectrum is illustrated in Fig. 1. Note 
that it is not necessary for the spectrum to have a (roton- 
like) minimum in order for there to be a critical velocity 
at non-zero momentum. We refer to excitations of mo- 
mentum near p c as levons. 

In the next section we lay out the order parameter de- 
scribing flow faster than the critical velocity, and calcu- 
late the corresponding energy of the system. As we show 
the order parameter develops an instability at the critical 
velocity towards formation of additional components at 
the critical momentum. In Sec. Ill we calculate, for flow 
velocities above the critical velocity for repulsive inter- 
excitation interactions and small amplitude modulation 
of the order parameter, the equilibrium order parameter, 
the spatially varying density induced by the additional 
components in the order parameter, the supercurrcnt, 
and the decrease in the superfluid mass density. We also 
qualitatively discuss the larger amplitude regime. In the 
final section we relate the present discussion to solitons 
in Bose-Einstein condensates, to vortices below the crit- 
ical angular velocity for vortex formation, and to super- 
critical currents in superconductors, as well as lay out a 
few possible ways that one could realize non-zero Landau 
critical velocities in cold atomic gases. 



II. ORDER PARAMETER WITH FLOW 

To examine the stability of the system we explore an 
order parameter in the lab frame, 



ip(x) 



a iqx 



(6) 



corresponding to a reduced number of particles in the 
initial condensate of momentum q, and with condensate 
components of momentum q±p c . The instability occurs, 
as we anticipate, at p c > 0, corresponding to an excita- 
tion moving in the negative x direction. We work at fixed 
average density, n. Such an order parameter, for small hi 
and V describes the usual Bogoliubov excitations of the 
condensate. While U and V can be generally complex, 
the most energetically favorable situation corresponds to 
U and V real, as will be apparent once we calculate the 
energy corresponding to the state Without loss of 
generality, we take U > 0, and U 2 > V 2 . 
The number density is 

n(x) = \tf;(x)\ 2 =n +U 2 + V 2 

+2^/n^{U — V) cosp c x — 2UV cos2p c x, (7) 

and the spatial average of the density is 



n = n +U 2 + V 2 . 



(8) 



Thus U 2 + V 2 is bounded above by fi. In the follow- 
ing we drop the bar over the n to simplify the notation. 
Similarly, the average kinetic energy density is 



E ke = (n q 2 +U 2 (q- Pc ) 2 +V 2 (q + p c ) 2 ) 
2m v ' 



2m 



(nq 2 + {U 2 + V 2 )p 2 c - 2{U 2 - V 2 )qp c ) .(9) 



The average interaction energy density is 



1 



l\\2 



Emt = 20 J d A Td A T'9{\T-r'\)\i>{x)\^{x')\ 

= ±g n 2 + gi n (U - V) 2 + U 2 V 2 g 2 , (10) 

where f2 is the system volume [l7| ■ We note that only the 
(U — V) 2 term in the energy distinguishes the relative sign 
of U and V, and this term, for given U and |V| favors V 
having the same sign as IA. In addition the mass current, 



mj = qn + (q - p c )U 2 + (q + p c )V 2 
= qn-{U 2 -V 2 )p c , 

is reduced from qn for p c > 0. 

To analyze the stability of the system we write 



(11) 



U = C cosh(?7/2), V = C sinh(r?/2), (12) 

where £ > 0, and < r\ < oo. Then the average number 
density is 



n = no + C 2 coshr;, 



(13) 



which implies the bound, C 2 cosh 77 < n. In terms of £ 
and 77 the average energy density is 

E = E ke + E mt = \ gm 2 + |-n + E', (14) 
2 2m 
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where 



Ef = C 2 (^cosh ?7 -M 
2m 



gine 



+C t.92 sinh r\ — g\e ''cosh 77 



(15) 



We suppress the two constant terms in Eq. (|14[) which 
play no role, and work with the energy density E' . 

Let us first ignore the £ 4 terms; minimizing the first 
line in (fT5j) with respect to 77 at fixed £, we find [l8j : 



{u-vf 
C 2 



2me c 



where 



2m V 2m 



2g\n 



1/2 



(16) 



(17) 



is the usual Bogoliubov expression for an excitation of 
momentum p c . For g\ > 0, 77 is positive. 
With this solution, 



E' = Q 2 



qpc 

m 



(18) 



We see explicitly that for q < mv c , the Landau critical 
momentum, the original supcrfluid flow is stable against 
small excitation as expected, but for q > mv c the system 
is unstable against developing a non-uniform condensate 
with C 7^ 0- The instability sets in, as noted, for p c > 0, 
corresponding to an excitation moving in the negative x 
direction. For ( < 1, the system energy E' is positive 
for all q < mv c . For q > mv c , however, the energy as 
a function of £ decreases steadily from zero indicating a 
transition to a state of nonzero £. 



III. ABOVE THE CRITICAL VELOCITY 

We now examine the state of the system for flow veloc- 
ity slightly above the critical velocity v c . While dE/dC, 2 
is initially negative, the terms of order £ 4 can, depending 
on parameters, locally stabilize the system at small ( 2 . 
To order £ 4 , we may take 77 to be given by Eq. ([TB]l. since 
the £ 2 corrections to 7/ do not contribute to E' at the 
minimum with respect to 77. The terms to order £ 4 in E' 
become 



E' 



where 7, given by 



C 2 s 



qpc 
m 



2 s 



9i 

2e 2 



5i.92" 



m V m 



(19) 



(20) 



corresponds, in terms of Pitaevskii's description of He-II 
above the Landau critical velocity Q , to the roton-roton 



interaction. Note that for small 52, 7 is negative, while 
for 



92 



> 



Pr 



- + 2<?in , 



gi m(gin) 2 \m 
7 > 0, and the system has a local minimum at 



c = (v-v c y- 



7 



(21) 



(22) 



This minimum corresponds to a nonzero condensate am- 
plitude at momentum p c . The excitations near p c form 
a levon condensate, or in the context of 4 He, a "Bose 
condensation" of rotons [f|. 

The variation of the energy with £ 2 about this mini- 
mum, which vanishes by construction, defines the renor- 
malized energy of an excitation of momentum p c in the 
lab frame, 



dE' 



+ C 2 7 



qp c 

m 



0. 



(23) 



This structure agrees with Pitaevskii's picture that the 
roton spectrum in the frame moving with the supcrfluid 
is shifted by the presence of critical excitations, e c — ¥ 
e c + £ 2 7 = e C) r, and the renormalized excitation energy 
in the lab frame, £ c ,r — IPc/m, vanishes. 

Furthermore the density has a "layered structure," 
Eq. ([7]). given to first order in ( by 



n(x) = n + 2vm~oC cosh 77 cos(p c x). 



(24) 



This structure, first predicted by Pitaevskii, has been 
seen in a density-functional simulation of flow in 4 He, 
with amplitude ~ (v — v^ 1 / 2 , corresponding to Eq. (|22|) 
Q. The total mass current in this state with nonzero £ 
is, from Eq. (|TT|) . 



mj = dE/dq — qn — p c Q 2 



(25) 



Thus slightly above the Landau critical velocity, the su- 
pcrfluid mass density is, 



mj p" c I v 

~ 77177, 1 

q 7 \Vc 



(26) 



Superfluidity does not cease at the Landau critical veloc- 
ity. Rather, as we see in this model calculation, as the 
Landau critical velocity is exceeded the supcrfluid mass 
density begins to fall below the total mass density, here 
linearly with the flow velocity The reduction in p s can be 
understood as a consequence of the lack of translational 
invariance of the system above v c 0]. Furthermore, the 
normal mass density is 



V c 7 \Vc 



(27) 



Quite generally, the supcrfluid mass density need not fall 
abruptly to zero at the critical velocity 
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The superfluid above the critical velocity has a richer 
mode structure than below. Ordinary zero sound - 
the conventional Bogoliubov excitations - should not be 
modified significantly for flows only slightly above v c . On 
the other hand, the presence of the normal component 
(the levon condensate) will lead to a additional sound-like 
mode, due to variations of the levon density and mean 
momentum. Noting that the calculations leading to Eq. 
(23) do not depend on q c being the critical momentum, 
we infer that the energy of an excitation of momentum p 
is given by 

w° = e v + C 2 7 - QP/m. (28) 

More generally, levons with momentum p are mixed with 
those with momentum 2p c —p because a pair of such exci- 
tations can be annihilated, thereby producing two levons 
with momentum p c . This process may also be considered 
as an Umklapp process resulting from the periodic lat- 
tice structure. As Pitaevskii showed Q for excitations of 
the "roton condensate" in superfluid 4 He, the excitation 
energies of the modes with momentum near the levon 
momentum for arbitrary direction of p are given more 
generally by 

up = [2C 2 7 (e P - q ■ p/m) + (e„ - q • P /mf)] 1/2 , (29) 

which reduces to the above result (f2"5|) for small positive 
v — v c , and p along q. The levon energy in the lab frame 
is 

^r.( V - Vc ) Pc (^- + 9l) (30) 

for small p - p c , where Q\\ = p x - p c , Qj_ = (p v ,p z ) 
and the /x's are effective masses. The structure of the 
levons is simply that of phonons of the layered conden- 
sate. Because the levon energy depends on the local par- 
ticle density and superfluid velocity, these modes will be 
hybridized with ordinary Bogoliubov zero sound excita- 
tions, thereby inducing an attractive interaction between 
levons. More detailed calculations of mode frequencies 
can be made using the methods employed for a Bosc- 
Einstein condensate in an optical lattice [H, , the only 
difference here being that the spatial modulation of the 
density is due to the flow, not an external potential. 

Let us look at the energy landscape more globally. For 
v slightly above v c and 7 > the system has a stable local 
minimum at small nonzero U and V. When 7 < 0, there 
is no local minimum at small £ 2 ; were the term of order 
£ 6 in the free energy positive, the system would undergo 
a discontinuous transition, leading to a discontinuity in 
p s . However, to describe the most favorable state of the 
system for 7 < requires that one goes beyond using 
the simple trial wave function, Eq. ([5]). Even for p c /2 < 
q < mv c , the energy E' initially rises with increasing £, 
indicating stability of the flow under small perturbations, 
but falls through zero at a nonzero £ determined from[l9j 
Eq. (fT5|) . On the other hand, for q < p c /2 the energy E' 
is positive for all nonzero £, indicating absolute stability 
of the flow. 



IV. DISCUSSION 

The wave function ([5]) for nonzero U and V, the low- 
est energy state at velocities greater than v c , may be 
regarded as a periodic array of solitons. Such arrays 
were first considered for the Gross-Pitaevskii equation by 
Tsuzuki and they have been studied for the case of a 
periodic background potential by Machholm et al. [Ioj |. 
who showed that they correspond to states on the upper 
branch of the "swallow tail" in the spectrum. In these 
works the solitons were found to have an energy greater 
than that of the uniform state moving at constant veloc- 
ity. The new feature of the situation considered in the 
present paper is that for velocities greater than w c , the 
periodic array of solitons has an energy less than that of 
the uniform state. This is a consequence of the nonzero 
range of the potential, which here makes density fluc- 
tuations with wavevectors close to p c less energetically 
costly than long-wavelength density fluctuations having 
the same amplitude. 

An approach similar to the present can be used to de- 
scribe instabilities of vortices in Bose condensates. Refer- 
ence [l2| considered a single vortex in a trap as a function 
of the rotation frequency, f2. For example, for il less than 
the critical frequency H c at which the system can first 
support a vortex at the center, the vortex is not stable 
against fluctuations, but has a negative energy anoma- 
lous mode [IH . The situation is analogous to that in the 
present paper above the critical velocity. The condensate 
wave function (cf. Eq. tp ~ y/me?^ + u + ve 2l< ^ , where 
4> is the azimuthal angle, then has lower energy, and de- 
scribes two off-center vortices in the cloud asymmetric 
about the origin. 

It is interesting to compare the present situation for 
bosons with the analogous phenomenon in superconduc- 
tors. As shown by Rogers [l4[ (quoted in [3]), when 
the superfluid velocity exceeds the Landau critical veloc- 
ity, v c = A/pp, where A is the gap and pf the Fermi 
momentum, the system spontaneously generates quasi- 
particles. At zero temperature each quasiparticle state 
of negative energy becomes occupied; the Pauli principle 
here plays the role of stopping the process from running 
away. Although the current continues to increase ini- 
tially with increasing flow velocity above v c , it quickly 
goes though a maximum and then goes to zero at veloc- 
ity (e/2)v c , where e = 2.718. . . is Napier's constant. In 
contrast to the Bose case, the excitations are fermions, 
so they cannot all be in the same momentum state. As 
a consequence, there is no density modulation for flow 
velocities greater than v c . 

In this article we have only studied states in which 
the amplitude of the spatial non-uniformity is small: in 
the future it would be interesting to extend this work 
to larger amplitude disturbances, including the possibil- 
ity that the order parameter could develop nodes, and 
even vortex lines could be formed. Numerical simula- 
tions in this non-linear regime would be very helpful. A 
further problem is to determine the state of the system 
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at finite temperature, including fluctuations of the con- 
densate structure from the levons. 

Finally, let us ask whether one can realize the present 
model experimentally. A first way could be to take ad- 
vantage of enhancement of the effective range near a 
Feshbach resonance. Essentially, for scattering of two 
particles at relative momentum p with total momentum 
0, the contribution of a resonance to the effective inter- 
action, g(p), is 



~l(p)r 



AE + p 2 /4m' 



(31) 



where AE is the energy difference of the particles of 
zero relative momentum in the entrance and intermediate 
channels, and \M\ 2 is the square of the matrix element 
between the two channels. Thus 



dg{p) 

l d(p 2 /2m) 



|M| 5 



2 (AE+p 2 /4m) 



(32) 



which can reach the needed value —1/2 sufficiently close 
to the Feshbach resonance at small p. A second possibil- 
ity could be to use the downward bending of the quasi- 
particle energy e{p) with increasing p, arising from short 
range correlations. To the extent that the f-sum rule 
is exhausted by a single quasiparticle, the quasiparticle 
energy is given by p 2 /2mS(p), where S(p) is the static 



structure factor [16| ; quite generally short range correla- 
tions tend to enhance S(p) at shorter wavelengths - in 
4 He leading to the roton dip, as Feynman first argued 
- and could produce the needed softening of the spec- 
trum. A third possibility could be to use atoms with 
dipole-dipole interactions in quasi-one dimensional sys- 
tems. We leave the elucidation of these methods to a 
future study. 

As we have shown in the calculations in this paper 
for weakly interacting condensed bosons, a uniform flow 
at velocity exceeding the Landau critical velocity does 
not necessarily destroy the superfluidity For positive 
interactions among the critical excitations, the system 
becomes stabilized with a reduced superfluid mass den- 
sity. A small amplitude analysis for negative interactions 
among the quasiparticles is however inadequate, and elu- 
cidating the possible states to which the system might go 
requires further work. 



Acknowledgements 

Author GB is grateful for the hospitality of the Aspen 
Center for Physics, where part of this work was carried 
out. We thank Tony Leggett and Sohcil Baharian for 
critical comments. This research was supported in part 
by NSF Grants PHY07-01611 and PHY09-69790. 



[1] L. D. Landau, J. Phys. USSR 5, 71 (1941). 

[2] L.P. Pitaevskii, J. Low Temp. Phys. 87, 127 (1992). 

[3] M. Greiner, O. Mandel, T. Esslinger, T. W. Hansch, and 

I. Bloch, Nature 415, 39 (2002). 
[4] L. P. Pitaevskii, Pi'sma Zh. Eksp. Theo. Fiz. 39, 423 

(1984) [Engl. Transl.: JETP Lett. 39, 511 (1984)]. 
[5] C. Raman, M. Kohl, R. Onofrio, D. S. Durfee, C. E. 

Kuklewicz, Z. Hadzibabic, and W. Ketterle, Phys. Rev. 

Lett/ 83, 2502 (1999). 
[6] S. V. Iordanskii and L. P. Pitaevskii, Usp. Fiz. Nauk 

131, 293 (1980) [Engl. Transl: Sov. Phys. Usp. 23, 317 

(1980)]. 

[7] A. Paramekanti, N. Trivedi, and M. Randeria, Phys. Rev. 

B 57, 11 639 (1998). 
[8] F. Ancilotto, F. Dalfovo, L. P. Pitaevskii, and F. Toigo, 

Phys. Rev. B 71, 104530 (2005). 
[9] B. Wu and Q. Niu, Phys. Rev. A 64, 061603 (2001). 
[10] M. Machholm, C. J. Pethick, and H. Smith, Phys. Rev. 

A 67, 053613 (2003). 
[11] T. Tsuzuki, J. Low Temp. Phys. 4, 441 (1971). 
[12] S. Baharian and G. Baym, Phys. Rev. A 82, 063606 
(2010). 

[13] D. S. Rokhsar, Phys. Rev. Lett. 79, 2164 (1997); M. Linn 
and A. L. Fetter, Phys. Rev. A 60, 4910 (1999). 



[14] K.T. Rogers, Ph.D. thesis, University of Illinois 1960. 
[15] J. Bardeen, Rev. Mod. Phys. 34, 667 (1961), Appendix 
B. 

[16] R. P. Feynman, Phys. Rev. 94, 262 (1954). 

[17] One cannot gain energy by allowing U and V to be com- 
plex: U — > Ue l8u , V — > Ve lS " . Since the phase difference 
of U and V can be absorbed in the choice of origin of x, 
we need only consider the more general order parameter 
ip(x) = e iqx + e tB (Ue~ ipcX - Ve 1 ^)}. The net ef- 

fect on the energy is to replace the (U — V) 2 in Eq. (|15p 
by U 2 + V 2 — 2UV cos 9. This final term is most negative 
however for = 0. Thus allowing phase variations of U 
and V cannot decrease the energy. 

[18] The quantity e _ '' //2 corresponds in Pitaevskii's calcula- 
tion [J] to the modulus of the matrix element of the den- 
sity operator between the uniform ground state and the 
state with one critical excitation. 

[19] The state of least energy for q > p c /2 is always at r\ = 
and no = 0, corresponding to a condensate of momen- 
tum q — p c , but this state is not germane to the present 
discussion. This result is independent of the interactions; 
inclusion of the gi term in the energy cannot provide 
absolute stabilization of the system at nonzero 



